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Abstract 

Analytical expressions of the one- and two-body terms in the cluster expansion of the 
charge form factors and densities of the s-p and s-d shell nuclei with N=Z are derived. 
They depend on the harmonic oscillator parameter b and the parameter j3 which origi- 
nates from the Jastrow correlation function. These expressions are used for the systematic 
study of the effect of short range correlations on the form factors and densities and of the 
mass dependence of the parameters b and (5. These parameters have been determined by 
fit to the experimental charge form factors. The inclusion of the correlations reproduces 
the experimental charge form factors at the high momentum transfers (q > 2 fm _1 ). It 
is found that while the parameter f3 is almost constant for the closed shell nuclei, 4 He, 
16 and 40 Ca, its values are larger (less correlated systems) for the open shell nuclei, 
indicating a shell effect in the closed shell nuclei. 

PACS numbers: 21.10.Ft, 25.30.Bf, 21.45. +v, 21.60.Cs 



1 INTRODUCTION 

The calculation of the charge form factors, F c h(q), and density distributions, p c h(r) of nuclei, 
is a challenging and appealing problem A possibility to face this problem is by means 
of an independent particle model. This approach, which is particularly attractive because of 
its simplicity, fails to reproduce the high momentum transfer data from electron scattering in 
nuclei |2], |, f|, |^, ^ [7|, [8], |9|, |H], |TTJ. For this reason a modification of the single particle (SP) 



potentials has to be suitably made. 

In fact a short range repulsion in this potential seems advisable for light nuclei [E2|. For 
example, with an harmonic oscillator (HO) potential having in addition an infinite soft core, 
the F ch (q) of A He can be well reproduced, but for the heavier nuclei, such as 12 C and 16 0, state 
dependent potentials seem necessary and even then the fit is not so good for higher q-values 



1 



12] . Another way is the introduction of a modified shell model with fractional occupation 



numbers for the various states [O, O, KB, [16j 



An approach, which is rather similar, is the introduction of the short range correlations 
(SRC) in the Slater determinant. Many attempts have been made in this direction, concerning 
mainly light closed shell nuclei in the framework of the Born-approximation [f|, |5|, || ^, |8|, [|, 



Czyz and Lesniak [[| first showed that the diffraction character of the 4 He form factor can be 
qualitatively explained by means of Jastrow-type [[2(J correlations. Khana ||, using Iwamoto- 
Yamada cluster expansion and retaining only one- and two-body cluster terms, showed that the 
inclusion of short range nucleon-nucleon correlations provide an adequate description for the 
known data on the elastic scattering of electrons by 40 Ca and make predictions for the behavior 
of the cross section at large momentum transfers. Ciofi Degli Atti, using the "single pair 
approximation" ||] and the Iwamoto-Yamada cluster expansion p] in s-p shell nuclei, showed 
that the elastic electron scattering at high momentum transfers seems to give a strong indication 
of the presence of SRC in nuclei. Bohigas and Stringari || and Dal Ri et al || evaluated the 
effect of SRC on the one- and two-body densities by developing a low order approximation 
(LOA) in the framework of the Jastrow formalism. They showed that one-body quantities like 
the form factor provide an adequate test for the presence of SRC in nuclei, which indicates that 
the independent-particle wave functions cannot reproduce simultaneously the form factor and 
the momentum distribution of a correlated system. Stoitsov et al [TTJ generalized the model 
of Jastrow correlations, suggested by Bohigas and Stringari [|J within the LOA of Ref . , to 
heavier nuclei like 16 0, 36 Ar and Ca reproducing very well the experimental data. 

In the above approaches, different types of expansions were used. The expansions were 
connected with the number of the simultaneously correlated nucleons, and the order of the 
Jastrow correlation function, /(ry), which were retained in the cluster expansion. Usually 
they are truncated up to the two-body terms which give significant contribution to various 
expectation values. 

In a series of papers [|17|, [Tj|, [OJ an expression of the elastic charge form factor, F ch (q), 
truncated at the two-body term, was derived using the factor cluster expansion of Clark et al 



21] , WA -, |23| . This expression, which is a sum of one-body and two-body terms, and depends on 
the HO parameter and the correlation parameter through a Jastrow type correlation function, 
is used for the calculation of F c h(q) of 4 He, 16 and 40 Ca and in an approximate way for the 
open s-p and s-d shell nuclei. 

The motivation of the present work is the systematic study of the effect of SRC on the s-p 



and s-d shell nuclei by completely avoiding the approximation made in earlier work |T7], [TR IIU 



for open shell nuclei, that is the expansion of the two-body terms in powers of the correlation 
parameter, in which only the leading terms had been retained. General expressions for the 
F c h(q) and p c h( r ) were found using the factor cluster expansion of Clark et al and Jastrow 
correlation functions which introduce SRC. These expressions are functionals of the SP wave 
functions and not of the wave functions of the relative motion of two nucleons as was the case 
in many previous works 0, [K], 0. Because of that, it is easy to extrapolate them to the case 
of open shell nuclei and use them either for analytical calculations with HO wave functions or 
for numerical calculations when more realistic SP wave functions are used. An advantage of 
the present method is that the mass dependence of the HO parameter b (with the presence 
of correlations) and the correlation parameter f3 can be studied. These parameters have been 
determined, for the various s-p and s-d shell nuclei by fit of the theoretical F c h(q) to the 
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experimental ones. It is found that while the parameter (3 is almost constant for the closed 
shell nuclei, 4 He, 16 O and 40 Ca, it takes larger values (less correlated systems) in the open shell 
nuclei, indicating a shell effect for the closed shells. The present method has also been used 
with wave functions of Skyrme-type interactions. In this case there is only one free parameter, 
the correlation parameter (3 and the two body-term is a small correction to the one body-term 
of the density. 

The paper is organized as follows. In Section II, the general expressions of the correlated 
form factors and density distributions are derived using a Jastrow correlation function. In 
section III the analytical expressions of the above quantities for the s-p and s-d shell nuclei, in 
the case of the HO orbitals, are given. Numerical results are reported and discussed in section 
IV, while an outline of the present work is given in section V. 



2 CORRELATED DENSITY DISTRIBUTIONS AND 
FORM FACTORS 

If we denote the model operator, which introduces SRC, by JF, an eigenstate <3> of the model 
system corresponds to an eigenstate 

* = .F$ (1) 

of the true system. 

Several restrictions can be made on the model operator JF, as for example, that it depends 
on (the spins, isospins and) relative co-ordinates and momenta of the particles in the system, it 



is a scalar with respect to rotations e.t.c. [p^| . Further, it is required that T is translationally 
invariant and symmetrical in its argument \ - ■■%■■■ A and possesses the cluster property. That 
is if any subset, %\ • • • i p , of the particles is removed far from the rest, i p+ \ • • -i^, T decomposes 
into a product of two factors, Ti\ ■ ■ ■ A) = T{%\ ■ ■ ■ i p ) JF(i p+1 • • • i A ) pjfl . In the present work 
T is taken to be of the Jastrow-type 

F = II finj) , (2) 

i<j 

where ffaj) is the state- independent correlation function of the form: 

/(r JJ ) = l-exp[-/3(r,-r J ) 2 ] . (3) 

The correlation function /(r^) goes to one for large values of =\ r 4 — | and it goes 
to zero for — > 0. It is obvious that the effect of SRC, introduced by the function /(r^-), 
becomes large when the SRC parameter (3 becomes small and vice versa. 

The charge form factor of a nucleus, in Born-approximation, can be written 

F ch (q) = f P (q) f DF (q) fcAi(q) F p (q) , (4) 

where f p {q) and fDF{q) are the correction for the finite proton size and the Darwin-Foldy rela- 
tivistic correction, respectively p5| , fc'Ai{q) is the Tassie-Barker [26[] center-of-mass correction 



and F p (q) the point form factor of the nucleus which is the expectation value of the one-body 
operator, 

A A 

°g = H°<?(A> = H ex P[«q r i] • ( 5 ) 



i=l i=l 
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That is, 



F P (q) = [ ^ ; = N(*\O q m = N(O q ) 



(6) 



where N = (^f(ri, r 2 , ■ ■ • , r^)|\l/(r 1 , r 2 , • ■ ■ , r^)) 1 is the normalization factor which is deter- 
mined so that F ch (0) = F p (0) = 1 or An / °° p(r)r 2 dr = 1. 
The point density distribution has the form 



P P ( r ) 



iV(tf|O r |tf) = JV(O r ) 



where 



A A 

O r = J2o r (i) = Xl'K 1 * -1 '*) 

i=l i=l 



(7) 



(8) 



2.1 EXPRESSIONS OF THE CORRELATED DENSITY DISTRI- 
BUTIONS AND FORM FACTORS 

In order to evaluate the point density distribution, p p {r), we consider, first, the generalized 
normalization integral, 

J(a) = (*|exp[a/(0)O r ]|*), (9) 
corresponding to the operator O r , from which we have, 



(Or) 



din I (a) 

da 



(10) 



a=0 



For the cluster analysis of equation (|T0|), following the factor cluster expansion of Ristig and 
Clark |2l], ^2], p3fl , we consider the sum-product integrals, hj((x), • • •, for the subsystems 

of the A-nucleon system and a factor cluster decomposition of these integrals. The expectation 
value of the density distribution operator is written in the form, 



(O r ) = (O r ) 1 + (O r ) 2 + --- + (O r ) A 



where 



(Orh 



E 

i=l 



da\Ii{a) 
da 



I ^i f o,(l)^i I i) 



(12) 



a=0 



i=l 



(°rh = E^[ ln ^'(«)- ln ^(«)- ln ^(«)Uo 

i<j OOL 

A A 

= I ^ 2 (or(l) + o r {2))T l2 | ij) a - £[<i I o r (l) I i) + (j I 0,(2) | j)} , (13) 

and so on. T\ is chosen to be the identity operator. 

The cluster expansion establishes the separation of one-body, two-body, • • •, A-body cor- 
relation effects on the density. Three- and many-body terms will be neglected in the present 
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analysis. Thus, in the two-body approximation, writing also the two sums of equation ( |13l) in 
a different form, p p (r) is written as follows, 



p p {r) = N(O r ) « iV[(O r )i + (O r ) 22 - (O r ) 21 ] 



where 



(O^ = J2(i\0 r (l)\i) , 

1=1 

A 

(O r } 2 2 = 2^{ij I Fl 2 o r {l)T 12 | ij) a 



i<j 
A 



(O r ) 21 = 2^(ij|o r (l)|ii) a . 

If the two-body operator is taken to be the correlation function given by (|3]), then 

•7^12 = 1 - 2(/(n, r 2 , /?) + ff (n, r 2 , 2(3) , 

where 

g>(ri, r 2 , z) = exp[— zrf\ exp[— zr%] exp[2zrir 2 coso;^] , z = (3,2(3 , 
and the term (O r ) 22 is written, 

(O r ) 22 = (O r ) 21 - 20 22 (r,/3) + 22 (r,2/3) , 

where 

A 

C 22 (r,2;) = 2^2(ij | o r {l)g{r x ,r 2 , z) | y) • 



So, p p (r) takes the form, 



p p (r) « iV[(O r )i - 20 22 (r,(3) + 22 (r,2(3)} 



(14) 

(15) 
(16) 
(17) 

(18) 
(19) 

(20) 

(21) 

(22) 



The terms (O r )i and 22 (r, z) and the density p p (r) can be expressed also in the convenient 
form: 

(O r ) l = p SD (r) , (23) 



22 (r, z) 



J g(r,r 2 , z)[p SD (r)psD(r 2 ) - p| D (r, r 2 )]dr 2 



p p (r) » JV 



Psr)(r) + / [g(r, r 2 , 2/3) - 2g(r, r 2 , /?)] P5s(r)p5Z)(r 2 ) - p 2 SD (r, r 2 



(24) 

dr 2 , (25) 



where Psd{ t i, t 2) is the uncorrelated density matrix associated with the Slater determinant, 



i=l 



The diagonal elements of this gives the one body density distribution, 

Pfifu(ri) = Psr)(ri,ri) . 



(26) 



(27) 



It should be noted that, a similar expression for p p (r), given by equation (p5|) , was derived 
by Gaudin et al. in the framework of LOA. In LOA the Jastrow wave function, of the 
nucleus was expanded in terms of the functions: g = Pir^) — 1 and h = /(r^) — 1 and was 
truncated up to the second order of h and the first order of g. This expansion contains one- and 
two-body terms and a part of the three-body term which was chosen so that the normalization 
of the wave function was preserved. Expression (^) of the present work has only one- and 
two-body terms and the normalization of the wave function is preserved by the normalization 
factor N. 

In the above expression of p p (r), the one-body contribution to the density is well known 
and given by the equation, 



(O r )i = Aj2vm(2l + l)J-^(r)^(r) , 

nl 47F 



(28) 



where rj n i is the occupation probability of the state nl (0 or 1 in the case of closed shell nuclei) 
and <p n i{r) is the radial part of the SP wave function. 

An expression for the two-body term is usually found by making a transformation to the 
relative and the center-of-mass coordinates of the two interacting nucleons 0, [TO, O. This is 



because the Jastrow function /(r^) depends on the relative coordinates of the two nucleons. 
Here, the expression for the two-body term, that is of the term 22 (r,z), will be found by 
expanding the factor exp[2zrir 2 cosco^] in spherical harmonics [p7, 2£|. That is 



expftzr^ cos u 12 \ = 2vr £ l^(2zrir 2 )K fc * mfc (Oi)Yfc mfe (02 



km i. 



(29) 



where 



U k (2zr 1 r 2 ) 



exp(2^rir 2 coso;i2)i : fc(cosa;i2)d(cosu;i2) 



2i k (2zr 1 r 2 ) 



(30) 



ik{x) is the modified spherical Bessel function. 

Using the algebra of spherical harmonics, the term: 022(^5 z) takes the form: 



22 (r,z) = 4 ]T W% ,.(2i i + l)(2/ i + l) x 



k=0 



z = (3,2(3 



(31) 



where 



An-ihrnhkr \ 



^0n 1 / 1 ( r )0n 3 « 3 ( r ) eX P[- 2;r2 ] X 



K 2 h ( r2 ) ^ 4 h ( r 2) exp [- zr\ ] i k (2zrr 2 )r% dr 2 



(32) 



Thus the expression of the term 22 (r,z) depends on the SP wave functions and so it is 
suitable to be used either for analytical calculations with the HO potential or for numerical 
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calculations with more realistic SP potentials. Expressions ( p8|) and ([H]) were derived for the 
closed shell nuclei with N=Z, where r] n i is or 1. For the open shell nuclei (with N=Z) we 
use the same expressions, where now: < rj n i < 1. In this way the mass dependence of the 
correlation parameter /3 and the HO parameter b can be studied. 

Finally, using the known values of the Clebsch-Gordan coefficients and provided that 
4Ur) = Mr), i.e.: A£ft$*(r,z) = A^fifafoz), equation for the case of s-p and 

s-d shell nuclei, takes the form: 

22 (r, z) = 4 [3A°r (r, z)rjL + &<Z°(r, z) - 6<°P(r, z)) rfc + ^(r, z)rjl + 

(95^i°(r, z) - ^A™?(r, z) - ^A™£ (r, z)) ^ + 

(lJLCi°(r, *) + lJMgg8°(r, *) - GA^r, *)) ^ + 

(20AZf(r, z) + 20A™(r, *) - 10A™(r, *)) ^ + 

(6Mg°»(r, z) + 60^g°(r, z) - 12^^, *) - lMg^r, *)) ^ + 

(4C°(r, *) + 4AiS8»(r, *) - 2^S8°(r, *)) + 

(l2A^°(r, z) + 12AjgS°(r, *) - 6A™(r, z)) Vlp rj 2s + 

(20^Sg°(r, *) + 2QAiSg°(r, z) - 10A^ 2 (r, *)) ^] . (33) 

The point form factor F p (q) can be derived in two equivalent ways. The first one is to follow 
the same cluster expansion as in the case of the density distribution and the second one is to 
take the Fourier transform of the density distribution p p {r), 

F P {q) = J exp[iq r]p p (r)dr . (34) 
In both cases, the form factor takes the following form: 

F p (q) « N[(O q ) 1 - 20 22 (q, (3) + 22 (q, 2/3)] . (35) 
In the above expression, the one-body term is given by the equation 

(0,>! =±J2Vni(2l + 1) / Ki(r)Mr)3o(qr)r 2 dr , (36) 



while the two-body term 22 (q, z) is given by the right hand side of equations (plQ and 
by replacing the matrix elements A^£*^ fc (r, z) by A™l\l™ l £ k (q,z) given by the equation, 



(pl^ir^^n^ir^expl-zrljit^znr^rldri , (37) 



where jo(x) and ik{x) are the spherical Bessel and the modified spherical Bessel functions 
respectively. 
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2.2 1-s COUPLING SCHEME 



In the cases of mean field potentials with 1-s coupling, the wave function of the nlj state has 
the form: 

VWi(ri) = <Mn) ]T {Irmsm^m^YPWxT^Xrtl) , (38) 

mi i m s 

where <p n ij{ri) is the radial part of the wave function of the state nlj and x™ s (l), Xt"*(l) are 
the spin and the isospin wave functions respectively. 

Following the same procedure as in sub-section [2.1| , the same expression 
been found. In this case the terms (O r )i and 022{r,z) have the form: 

(O r ): =2Y,rhij{2j + l)^y(0<MO (39) 

nl 

and 



for p p (r) has 



<9 22 (r,z) = 2 £ (2j i + l)(2j j + l) x 

-1 



2 Anjljjjnikjik , 



fc=0 



z = 13,2/3 ,(40) 



where the matrix elements ^nlilf^hjt ( r > z ) are gi ven again by equation ( |32"D replacing the 
wave functions n i( r ) by the wave functions (f) n ij(r). Equation flip]) can also be written in a 
similar form with that of equation ( p3[ ) with the difference that there are now more terms. This 
expression goes to (p3|) if there is no 1-s coupling. 

In the evaluation of the density matrix Psd{^Xi ^2)-, which is necessary for the derivation of 
the expression ( (40|) of the term 022(f) 2), we made the following approximation: In the sum 
over the spin coordinates, only terms of pairs of particles having the same third spin component 
are taken into account. The contribution of the terms which contain pairs with opposite third 
spin component is small and is neglected [^. [30|| . In this scheme the form factor is calculated 
numerically by Fourier transform of p p (r) employing equation (|34T). 



3 ANALYTICAL EXPRESSIONS 

In the case of the HO wave functions, with radial part, 



Mr) = N nl b- 3/2 




cxp 



2b 2 



2n\ 



1/2 



nl 



T(n + l + 



(41) 



where the Associated Laguerre polynomial is defined as in Ref. |JT|, analytical expressions of 
the one-body term and of the matrix elements A"^"*j* fe (r, z) and A^^f (q, z) , defined by 
equation fl32f ) and fl3~7j) , can be found. From these expressions, the analytical expressions of the 
terms 22 (r, z) and 22 (g, z), defined by equation (|33|), can also be found. 

The expression of the one-body term of the density and form factor has the form: 



C e 



2 

E 

fc=0 



2k 



X 



r,q 



(42) 
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where for the (O r )i the variable £ and the coefficients C and C 2 k are: 

2 



Co 
C 2 

c 4 



c 



r 

2t7i s + 3r? 2s , 

o ( 2? 7ld + V2s) , 



(43) 



while for the (0 9 )i the corresponding quantities are: 



Co 
C 2 

C 4 



C 



2(Vn + V2s + 3i]i P + 5i] 



id) 



(44) 
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(2ry w + r]2s 



The analytical expression of the matrix element A^f?™^ (r, z), which are given by equation 
J), has the form: 



n ^ 

\i=l 

EE " 

w=0 s=0 

{\{k + i 



y 



I 16v^F6 3 

ir +s (n 2 + l 2 + l 
n 2 — w 



1 + 2j/ 2 



1+y 



n 4 + Z 4 + | 



wis! 



n 4 — s 



4 -k) + W + S)\ 
1 _|_ y\^(h+h+k+3)w+s \(l 2 +li-k)+w+s 



y2 e 



l + y 



(45) 



where £ = r/b and £/ = zb 2 (z = j3, 1(5). It is mentioned that the Clebsch-Gordan coefficients 
(Ii0lj0\k0) , which appear in equation (|31~|) (where k runs from to /j + Zj), are different from 
zero only when Zj + Zj + is an even integer and therefore l 2 + Z 4 — /c is an even integer. Because 
of that the lower index of the Laguerre polynomial of equation fl4~5|) is an integer. 



The above expression of A$£fc k (r,z) is of the form: /(£ 2 ) exp[-ig^ 2 ], w h ere /(£ 2 ) is 



a polynomial of £ 2 . The substitution of the expression of A^f^f^fr, z) to the expression of 
22 (r,z), which is given by equation leads to the analytical expression of the two-body 
term of the density. This expression, for the case of s-p and s-d shell nuclei, is again of the 
form ,/(£ 2 ) exp[— i^£ 2 ] where now /(£ 2 ) is a polynomial of the form: 



/(e a ) = £c»fo»*i*,v)£ 



2 k 



k=0 

The corresponding analytical expression of the matrix element A™^™*^ (q, z) 



(46) 
which are 



given by equation (|37|) , has the form: 



\i=l 



i + y ,2 

l + 2y . 



ni n 3 712 n,4 2 



\{h+h-k)+w+s 



EEEE E 

p=0v=0u)=0s=0 t=0 



/_]\p+v+uf+s+t 

p\v\w\s\tl 
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T^(l 2 + U + k + 3)+w + s ] J 


(h(h+h + k) +p + v + t 




r(t + fc + §) 





f 1 \ ll \\(h+H—l2—l4)+p+v—w — s i f -l i 

„2t+fc I 1 + ?/J 2 7-2 / _LZA<2 \ , , 7) 



(1 + 2y)!( /l+ ' 3+fc+3 ) +p+v+ * |('i+'3+fe)+p+v+t VI + 2y 

where £ = g&/2 and y = z& 2 (2; = (3, 2/5). The lower index of the Laguerre polynomial of 
equation fl4"?|) is an integer because Zi + / 3 + k is an even integer. 

The above expression of A^^f (q, z) is of the form: /(£ 2 ) exp[-i^£ 2 ], where /(£ 2 ) is 

a polynomial of £ 2 . The substitution of the expression of A^ff™* 1 ^^, z) to the expression of 
022(q,z), given by equation (|33|) , leads to the analytical expression of the two-body term of 
the form factor. This expression, for the case of s-p and s-d shell nuclei, is again of the form 
/(£ 2 ) exp[— y^£ 2 ] where now /(£ 2 ) is a polynomial similar to that of equation (|46D. 

From the analytical expression of the form factor, the analytical expression of the mean 
square radius of a nucleus, which is the coefficient of —q 2 /6, can be found. This expression is 
of the form: 

(r 2 ) = iVKr 2 )! - 2(r 2 (/% 2 )> 22 + (r 2 (2/% 2 )> 22 ] , (48) 



where 
and 



(r 2 }i = [6(77!, + V2s) + 30r7 lp + 70r] ld ]b 2 , (49) 
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(r\y))22 = b\l + 2y)- 13 / 2 £ a k y k , y = [3b 2 , 2(3b 2 . (50) 

fc=0 

The coefficients aj~ depend on the occupation probabilities of the various states. If we 
expand, the right hand side of equation (0), in powers of 1/y (y = fib 2 > 1) and keep powers 
of f3b 2 up to (/36 2 ) -3 / 2 , an approximate expression of the contribution of the two-body term to 
the radius takes the simple form: 

(r 2 (y)) 22 = Cb 2 ((3b 2 )- 3 / 2 , (51) 

where the coefficient C depends on the occupation probabilities of the various states. 

The analytical expressions of the form factor and the density, which were found previously, 
will be used in section IV for the fit of the theoretical charge form factors to the experimental 
ones and for the calculations of the charge density distributions for various N=Z (s-p and s-d 
shell) nuclei. 



4 RESULTS AND DISCUSSION 

The calculations of the charge form factors for various s-p and s-d shell nuclei, with N=Z, have 



been carried out on the basis of equations (|J) and (|35| ) and the analytical expressions of the one- 
and two-body terms which were given in section III. Two cases have been examined, named 
Case 1 and Case 2, which correspond to the analytical calculations with HO wave functions 
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without and with SRC respectively. In Case 1 there is one free parameter, the HO parameter 
b, while in Case 2 there are two free parameters, the parameter b and the correlation parameter 
p. The parameters, in both cases, have been determined, for each nucleus separately, by least 
squares fit to the experimental F c h(q). 

The best fit values of the parameters as well as of the values of x 2 , 

X 2 = ^J2m(the) - F t (exp))/AF t (exp)] 2 , (52) 

A i=l 

are displayed in Table I. In the same table the calculated root mean square (RMS) charge radii 
( r ch) 1//2 an d the contribution of the SRC to them, 

(r 2 h = (r 2 ch ) - (r^/A , (53) 

are displayed and compared with the corresponding experimental RMS radii. It is noted that 
(r 2 )2 is independent from the centre-of-mass correction and finite proton size. 

The experimental and the theoretical F c h(q), for the various cases, for the closed shell nuclei: 
4 He, 16 and 40 Ca are shown in Figure 1 while for the open shell nuclei: 12 C, 24 Mg, 28 Si and 
32 S, are shown in Figure 2. 

From the values of \ 2 ■> which have been found in Cases 1 and 2 (see Table I) and also from 
Figures 1 and 2, it can be seen that the inclusion of the correlations improves the fit of the form 
factor of all the nuclei we have examined. Almost all the diffraction minima which are known 
from the experimental data are reproduced in the correct place. There is a disagreement in 
the fit of the form factor of the open shell nuclei 24 Mg, 28 Si and 32 S for q « 3.5 fm -1 where 
it seems that there is a third diffraction minimum in the experimental data, which cannot be 
reproduced in both cases. 

It is seen from Table I that the parameter b has the same behavior as function of the mass 
number A in the HO and the correlated model, while the following inequality holds: 

b(HO) > b(SRC) . 

This is due to the fact that the introduction of SRC tends to increase the relative distance 
of the nucleons i.e. the size of the nucleus, while the parameter b, which is (on the average) 
proportional to the (experimentally fixed) radius of the nucleus, should become smaller. Such 
a behavior should be expected in view also of relations fl49|) and (|51~1) . 
It is also noted that the difference: 

Ab = b{HO) - b{SRC) , 

is almost constant for the open shell nuclei and it is larger for the closed shell nuclei 4 He, 16 
and 40 Ca. This can also be seen from figure 3 where the values of Ab versus the mass number 
A has been plotted. The behavior of Ab as function of A indicates that the SRC are stronger 
for the closed shell nuclei than in the open shell ones. 

In Figure 4 the values of the correlation parameter {3 versus the mass number A have been 
plotted. From this figure it is seen that the parameter (3 is almost constant for 4 He, 16 0, and 
40 Ca and takes larger values (less correlated systems) in the open shell nuclei. 

The behavior of the two parameters, b and /3, indicates that there should be a shell effect 
in the case of closed shell nuclei. That is, there is a shell effect not only on the values of the 
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harmonic oscillator spacing hu, as has been noted in Refs. [37, |38| but also on the values of 
the correlation parameter (3. 

In the above analysis, the nuclei 24 Mg, 28 Si and 32 S were treated as Id shell nuclei. We have 
also considered the Case 2* in which the occupation probability r]2 S of the nuclei 24 Mg, 28 Si and 
32 S is taken to be a free parameter besides the other two parameters b and f3. We found that 
while the \ 2 values become better, comparing to those of Case 2, the third diffraction minimum 
is not reproduced either and the behavior of the parameters b and (3 as functions of the mass 
number A remains the same. The results in this case are shown in Table I and in figure 2. 
The values of the occupation probability r]2 S of the above mentioned three nuclei are: 0.0355, 
0.0245 and 0.2945 respectively, while the corresponding values of rjid, which can be found from 
the values of r] 2s through the relation: 

Vld = [{Z - 8) - 277 2s ]/10 

are: 0.3929, 0.5951 and 0.7411 respectively. The partial occupancy of the state 2s for the 
nucleus 32 S has as result to increase the central part of the charge density significantly. 

In Figures 5 and 6 the charge densities, p c h{ r ) (normalized to Z), of the above mentioned 
nuclei for the various cases of Table I are shown and compared with the charge densities of Ref. 



[32j. In the same figures the contribution of the SRC to p c h(r), 

P2,ch(r) = Pch{r) - Pi,ch(r) , (54) 

is shown. The introduction of short range correlations has the feature of reducing the central 
part of the densities of the closed shell nuclei, while P2,ch{ r ) is small and characterized by 
oscillations in the case of open shell nuclei. 

From the determined mass dependence of the parameters b and (3, the values of these 
parameters for other s-p or s-d shell nuclei can found. In Figures Id and 5d the F c h(q) and 
Pch{ r ) of the nucleus 36 Ar, treated as an Id closed shell nucleus, are shown. As there are no 
experimental data for F c h(q) for high q values, the value of the parameter (3 is taken to be the 
mean value of the corresponding values of 16 and 40 Ca, that is : (3 3 q = 2.2937 fm -2 , while the 
parameter b is determined assuming that 

A6 36 = Ai> 40 , 

where 

Ab A = b A (HO) - b A (SRC) . 

Using the values of the parameters b^HO) = 1.9453 fm, b^SRC) = 1.8660 fm from Table I 
and choosing the parameter b 3e (HO) = 1.8800 fm in order to reproduce the experimental RMS 
charge radius of 36 Ar (< r 2 > l J x 2 p = 3.327 ± 15 fm |32|) the value b 36 (SRC) = 1.8007 fm was 
found. These values of /? 36 and b 36 (SRC) have been used for the calculations of the correlated 
F c h(q) and p c h{f) of 26 Ar, which are shown in figures Id and 5d, respectively. The calculated 
RMS charge radius, < r 2 > 1 / 2 = 3.3343 fm, which was found, is within the experimental error. 

The effect of SRC on the form factors has also been examined with various Skyrme-type 
wave functions. In this case, there is only one free parameter, the correlation parameter (3. The 



potential parameters have been adjusted in Refs. |39]. [10| in order to reproduce various physical 
quantities such as RMS radii, binding energies e.t.c. Thus, some effects of the SRC, we would 
like to study, have already been averaged out. Because of that we cannot study the effect of 
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the SRC on the parameters of the Skyrme interactions. This can be done if both the Skyrme 
parameters and the correlation parameter are readjusted by fit to various experimental data 
including the form factors and the fit is made for many nuclei at the same time. As this is out 
of the scope of the present work, we examined the form factors of 16 and 40 Ca with various 



Skyrme interactions, namely SKI to SK6 |39], f|(J, without SRC, named Case 3 and with SRC, 
named Case 4. 

In Case 3 and for 16 0, only SKI gives smaller value for x 2 compared with that of Case 1 
(HO without SRC). The inclusion of SRC (Case 4) to the Skyrme-type wave functions gives 
better \ 2 onr y f° r the SKI, but still this value is about 10% larger than in Case 2 (HO with 
SRC) and 1% smaller than in Case 3. For SK2 to SK6 the correlated parameter j3 goes to very 
large values (very small correlations) without improving the quality of the fit. For 40 Ca and for 
Case 3 all the Skyrme interactions, we have examined, give almost the same \ 2 value which is 
about 20% smaller than in Case 1 and 10% larger than that of Case 2, while the inclusion of 
the correlations improves the quality of the fit, for all the Skyrme interactions, not more than 
2%. See Table I and figure 1 for the results we have found with SKI. 

From the results with Skyrme-type wave functions we could conclude that even if a mean 
field is more realistic than the HO one, the inclusion of the SRC does not improve the fit, at 
least of the F ch (q), significantly. Significant improvement of the fit should be expected if the 
parameters of the mean field will be determined along with the correlation parameter. It should 
be noted also that, a good fit to the experimental F c h(q) of 6 Li, 12 C and 16 has been found by 
Ciofi Degli Atti et al , using Woods-Saxon wave functions, only if the correlation parameter 
and the depths and radii of the wells were free parameters . 



5 SUMMARY 

In the present work, general expressions for the correlated charge form factors and densities 
have been found using the factor cluster expansion of Clark et al. These expressions can be 
used, either for analytical calculations, with HO orbitals or for numerical calculations of F c h(q) 
and Pchir), with more realistic orbitals. 

The systematic calculations, with HO orbitals, of these quantities, for a number of N=Z, 
s-p and s-d shell nuclei, show that there is a shell effect on the values of the HO parameter b 
and on the correlation parameter f3. Regarding the parameter (3 it is almost constant for the 
closed shell nuclei while it takes larger values for the open shell nuclei. The mass dependence 
of these parameters indicates that the SRC are stronger for the closed shell nuclei than for the 
open shell ones. This dependence can also be used to find reasonable values of these parameters 
for nuclei for which there are no experimental data of the F c h(q). 

Numerical calculations with various Skyrme-type wave functions, taken from the literature 



39], |40[1 , indicate that even if a mean field is more realistic than the HO one, the introduction 
of the SRC does not improve the fit of the charge form factors significantly. The reason for this 
should be that the introduction of the correlations makes a change in the mean field in a way 
that there is a balance between the SRC and the mean field, while for a given mean field there 
is not this flexibility. A way to overcome this difficulty should be to readjust the parameters of 
the mean field and the parameter of the correlations. 



13 



ACKNOWLEDGMENTS 



The authors would like to thank Professor M.E. Grypeos and Dr. CP. Panos for useful com- 
ments on the manuscript. 

References 

[1] a. L.R.B. Elton, Nuclear Sizes (Clarendon Press, Oxford 1961). 

b . H. Uberall, Electron Scattering from Complex Nuclei (Academic Press, New York and 

London 1971). 

c. R.C. Barrett, and D.F. Jackson, Nuclear Sizes and Structure (Clarendon Press, Oxford 

1977). 

[2] W. Czyz, and L. Lesniak, Phys. Lett. B 25, 319 (1967). 
[3] F.C. Khana, Phys. Rev. Lett. 20, 871 (1968). 
[4] C. Ciofi degli Atti, Nucl. Phys. A129, 350 (1969). 
[5] C. Ciofi degli Atti, and N.M. Kabachnik, Phys. Rev. C 1, 809 (1970). 
[6] G. Ripka, and J. Gillespie, Phys. Rev. Let. 25, 1624 (1970). 
[7] M. Gaudin, J. Gillespie, and G. Ripka, Nucl. Phys. A176, 237 (1971). 
[8] O. Bohigas, and S. Stringari, Phys. Lett. B 95, 9 (1980). 
[9] M. Dal Ri, S. Stringari, and O. Bohigas, Nucl. Phys. A376, 81 (1982). 
[10] H.P Nassena, J. Phys. G 14, 927 (1981). 

[11] a. M.V. Stoitsov, A.N. Antonov, and S.S. Dimitrova, Z. Phys. A 345, 359 (1993). 
b. M.V. Stoitsov, A.N.Antonov, and S.S. Dimitrova, Phys. Rev. C 47, 2455 (1993). 

[12] a. M. Grypeos, and K. Ypsilantis, J. Phys. G 15, 1397 (1989). 

b. K. Ypsilantis, and M. Grypeos, 4th Hellenic Symposium on Nuclear Physics, Ioannina- 
Greece, (1993), 99 and J. Phys. G: 21, 1701 (1995). 

[13] B.A. Brown, S.E. Massen, and PE. Hodgson, J. Phys. G 5, 1655 (1979). 

[14] F. Malaguti, A. Uguzzoni, E. Verodini and P.E. Hodgson, Rivista Nuovo Cimento 5, 1 
(1982). 

[15] a. I.S. Gul'karov, Sov. J. Part. Nucl. 19, 149 (1988). 

b. I.S. Gul'karov, and B.P Nigam, Phys. Rev. C 52, 663 (1995). 

[16] T.S. Kosmas, and J.D. Vergados, Nucl. Phys. A536, 72 (1992). 

[17] S.E. Massen, H.P. Nassena, and CP Panos, J. Phys. G 14, 753 (1988). 

14 



[18] S.E. Massen, and CP. Panos , J. Phys. G 15, 311 (1989). 
[19] S.E. Massen, J. Phys. G 16, 1713 (1990). 
[20] R. Jastrow, Phys. Rev. 98, 1497 (1955). 

[21] J.W. Clark, and M. L. Ristig, Nuov. Cim. LXXA 3, 313 (1970). 

[22] M.L. Ristig, W.J. Ter Low, and J.W. Clark, Phys. Rev. C 3, 1504 (1971). 

[23] J.W. Clark, Prog. Part. Nucl. Phys. 2, 89 (1979). 

[24] D.M. Brink, and M.E. Grypeos, Nucl. Phys. A97, 81 (1967). 

[25] H. Chandra, and G. Sauer, Phys. Rev. C 13, 245 (1976). 

[26] L.J. Tassie, and F.C. Barker, Phys. Rev. Ill, 940 (1958). 

[27] R.R. Roy, and B.P. Nigam, Nuclear Physics, John Wiley & Sons, Inc 1967. 

[28] A. de-Shalit, and I. Talmi, Nuclear Shell Theory, Academic Press (1963). 

[29] F. Arias de Saavedra, G. Co', A. Fabrocini, and S. Fantoni, Nucl. Phys. A605, 359 (1996). 

[30] F. Arias de Saavedra, G. Co', and M.M. Renis Phys. Rev. C 55, 673 (1997). 

[31] I.S. Gradshteyn, and I.M. Ryzhik, Tables of Integrals, Series, and Products, Academic 
Press (1980). 

[32] H. De Vries, CW. De Jager, and C. De Vries, Atom. Data and Nucl. Data Tables, 36, 495 
(1987). 

[33] a. R. Frosch, Phys. Lett. B 37, 140 (1971). 

b. R.G. Arnold, B.T. Chertok, S. Rock, W.P Schutz, Z.M Szalata, D. Day, J.S. McCarthy, 
F. Martin, B.A. Mecking, I. Sick, and G. Tamas, Phys. Rev. Lett. 40, 1429 (1978). 

[34] I. Sick, and J.S. McCarthy, Nucl.Phys. A150, 631 (1970). 

[35] B.B. Sinha, G.A. Peterson, R.R. Whitney, I. Sick, and J.S. McCarthy, Phys. Rev. C 7, 
1930 (1973). 

[36] G.C. Li, M.R. Yearian, and I. Sick, Phys. Rev. C 9, 1861 (1974). 

[37] C.B. Daskaloyannis, M.E. Grypeos, C.G. Koutroulos, S.E. Massen, and D.S. Saloupis, 
Phys. Lett. B 121, 91 (1983). 

[38] a. G.A. Lalazissis, and CP Panos, Phys. Rev. C 51, 1297 (1995). 

b. G.A. Lalazissis, and CP. Panos, Intern. Journal of Mod. Phys. E 5, 664 (1996). 

[39] D. Vautherin, and D.M. Brink, Phys. Rev. C 5, 626 (1972). 

[40] M. Beiner, H. Flocard, N. Van Giai, and P. Quentin, Nucl. Phys. A238, 29 (1975). 



15 



Table 1: The values of the parameters b and /3, of the x 2 > and of the RMS charge radii {r^) 1 ^ 2 : 
contribution of the one body density (column HO), contribution of SRC (column SRC) and of 
the total RMS charge radii (column Total), for various s-p and s-d shell nuclei, determined 
by fit to the experimental F C h(q). Case 1 refers to the HO form factor, Case 2 when SRC are 
included. Case 2* is the 2 but with the occupation probability of the state 2s 

taken to be a free parameter. Cases 3 and 4 refer to SKI potential with and without SRC, 
respectively. The experimental RMS charge radii are from Ref. ||32j|. 
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Figure 1: The charge form factor of the nuclei: 4 He (a), 16 (b), 40 Ca (c) and 36 Ar (d) for 
various cases. The experimental points for 4 He, 16 0, and 40 Ca are from Refs. P5] 
respectively. 
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Figure 2: The charge form factor of the nuclei: 12 C (a), 24 Mg (b), 28 Si (c) and 32 S (d) for 
various cases. The case HO+SRC* corresponds to the case when the occupation probability 
rj2s is treating as free parameter. The experimental points of 12 C are from Ref. pi and for the 



other nuclei from Ref. [36|. 
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Figure 3: The difference Ab = b(HO) — b(SRC) versus the mass number A. b(HO) and b(SRC) 
are the HO parameter in Cases 1 (HO without SRC) and 2 (HO with SRC) respectively. 
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Figure 4: The correlation parameter (5 versus the mass number A in Case 2 (HO+SRC). 
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Figure 5: The charge density distribution p c / l (r) and the contribution P2,ch(j) of SRC to it of 
the nuclei: 4 He (a), 16 (b), 40 Ca (c) and 36 Ar (d) for various cases. The experimental points 
are from Ref. |32|| . 
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Figure 6: The charge density distribution p c h(j) and the contribution P2,ch(r) of SRC to it 
of the nuclei: 12 C (a), 24 Mg (b), 28 Si (c) and 32 S (d) for various cases. The case HO+SRC* 
corresponds to the case when the occupation probability r)2 S is treated as a free parameter. The 
experimental points are from Ref. |32 
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